Scale-free properties of weighted networks 
with connectivity-driven topology 



W. Jezewski 



1 



Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego 17, 

60-179 Poznan, Poland 



Abstract 

The rate equations are used to study the scale-free behavior of the weight distri- 
bution in evolving networks whose topology is determined only by degrees of pre- 
existing vertices. An analysis of these equations shows that the degree distribution 
and thereby the weight distribution remain unchanged when the probability rate of 
attaching new nodes is replaced with an unnormalized rate determined by the ratio 
of the degree of a randomly selected old node to the maximal node degree at the 
current stage of the network evolution. Such a modification of the attachment rule 
is argued to accelerate considerably numerical simulations of both unweighted and 
weighted networks belonging to the class of investigated evolving systems. It is also 
proved that the studied rate equations have a solution corresponding to the total 
weight (concentrated at individual vertices) distribution displaying the power-law 
behavior for asymptotically large weights. 
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1 Introduction 



Scale-free networks belong to an important group of complex evolving systems 
exhibiting a power- law of the degree (connectivity) distribution [1]. To explain 
the origin of the scaling form of the degree distribution, several scenarios of the 
growth and self-organization of networks have been proposed [2-11]. Probably 
the most natural mechanism responsible for scaling properties of time-evolving 
networks is the preferential attachment of a new vertex with the probability 
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being proportional to the degree of a randomly selected preexisting node to 
which a new one will be connected [2-4]. This simple self-organization mech- 
anism allows one to recover the complex topological structure of scale-free 
networks [1,3]. However, in many real network systems, links perform some 
specific functions, and the underlying binary networks can only be considered 
as skeletons determining topological properties of these systems. Thus, in real 
networks, links play a role of arteries through which some processes of different 
intensities are executed, or which carry loads of different magnitudes [12-17]. 
Functional properties of networks can be characterized by assigning loads or 
weights to particular links [18]. Then, the network evolution is in general de- 
scribed by both preferential node attaching and assigning weights to links. 
Essentially, processes of the preferential node attachment and the link loading 
are not autonomous. The preferential mechanism of the vertex attachment can 
be determined only by degrees of existing vertices [18], but can also be affected 
by weights, or can even be controlled entirely by total weights concentrated 
at preexisting vertices [15-17]. On the other hand, weights ascribed to newly 
created links can depend on degrees of old nodes to which new nodes are con- 
nected [18], or, in cases when internal (existing) links are loaded, the weight 
assigned to a given internal link can be expressed by degrees of both vertices 
joined by this link [12]. Models in which the attachment rule is entirely gov- 
erned by degrees of existing nodes appear to be appropriate for describing not 
only topological but also functional features of the WWW. In this network, the 
decision about linking a pair of nodes (webpages) depends on the popularity 
(measured by the node degree) of a given preexisting node to which an another 
node will be connected rather than on the weight concentrated at the existing 
node, i.e., the total size of documents sent from this node to others. Numerical 
simulations of networks with degree-driven topologies indicate that the total 
weight distribution can display in such systems a power-law form [18,19]. In 
turn, models whose topologies are entirely controlled by weights provide a 
representation of real systems in which a change of network traffic due to 
connecting a new node to an old node strongly disturbs functions performed 
by all links joining the old node with other nodes, or in which connections 
of new vertices influence functions of whole systems. Illustrative examples of 
such real networks are transportation systems [15-17]. It has been shown an- 
alytically that, in the case of linear dependence of the attachment probability 
on total weights concentrated at vertices, the distribution of total weights can 
also display a power-law form [15-17]. 

In this paper, undirected weighted networks with connectivity driven topolo- 
gies are considered. By analyzing the rate equations for the degree and weight 
distributions in such networks, it is shown that these distributions remain 
unchanged as the preferential mechanism of node attaching is appropriately 
modified by using an unnormalized attachment rule. This change of the ver- 
tex attaching rule is shown to be very useful for numerical simulations of 
some unweighted (binary) networks as well as for simulations of the consid- 
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ered weighted networks with degree-governed topological structures. The rate 
equations for the time evolution of the distribution of total weights concen- 
trated at vertices of a given connectivity is proved to have a solution of a 
scale-free form. 



2 Evolution of the degree distribution in binary networks 



In this section, the preferential attachment mechanism that a new node is 
linked to a randomly chosen preexisting node with the probability depending 
on the degree of the earlier node is revised. In general, this probability is 
assumed in the form [4,20] 

n^k) = Ay(k)/M^t) (1) 



with the connection kernel 



A 1 {k) = < 



(2) 

a k k ,7 = 1, 



and with the normalization factor 1/M y (t) determined by 

^max 

M 7 (f) = £ A,(j) Nfat) , (3) 

3=1 



where 7 > is the nonlinearity exponent, k denotes the vertex connectivity, 
kmax = ^max(7; t) is the maximal node degree, the amplitudes a k > 0, k = 
1,2,..., and iV fc (7;t) denotes the number of vertices of the degree k. Note 
that the power-law N k ~ k~ v emerges only if the kernel A^(k) is linear or 
asymptotically linear, i.e., when 7 = 1 and the amplitudes a k are all equal, or 
are unequal, but such that a k — > as k — > 00 [20]. 

It is proved below that the degree distribution remains unchanged when the 
probability attachment rate i7 7 (/c) is replaced with the rate 



afc&/aoo&max(l;i) , 7 = 1 . 



For cases of < 7 < 1, this modified attachment rate is unnormalized (thus, 
it does have the character of a probability rate), and is not proportional to 
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the rate II 1 (k) as, in these cases, time dependences of /c m ax(7;£) and M 1 {t) 
are different. If a k > (for some k < k max ), the rate TI-i{k) takes even 
values greater than one. To show the equivalence of both preferential linking 
procedures, one can apply the rate equations describing, for the attachment 
rate (1), the time evolution of the degree distribution [20] 



Solutions of these equations imply for large t that [20,21] 



M 7 (t) 




(6) 



where /i is a function of 7, such that, in cases of homogeneous kernels with 
uniform, unit amplitudes, ji G [1,2] for all 7 > (/i — 1 for 7 > 1), and 



&max(7; t) = < 



(1=7)1^ 

, . 1 



(In t) 



1-7 



, 0< 7 <! , 
, 7 = 1 , 
, 7 > 1 , 



(7) 



with v > 2, s being dependent on v as well as on the manner of the convergence 
of a*; to as k — > 00 (in the case of 7=1). Since, fc max (7;£) is a monotonic 
function of time for each 7, Eq. (5) can be transformed from t to fc max using 
the inverse function t = <7 7 (&; max ). Thus, the use of Eqs. (6) and (7) gives 

M~ f (t) = kl ax g'Jk max ) , (8) 



where, g' denotes the derivative of g 1 with respect to k max . Next, inserting (8) 
into Eq. (5), replacing t by A; max , and multiplying these equations by ^(fc max ), 
one obtains 



-N k (r, k max ) = 



rlk 

"■"■max 

— |— [A 7 (k - 1) ^-1(7; /c max ) - i 7 (/c) iV fe (7; /c max ) ] + S kA ^(A; max ) (9) 

fvmax 

with N k (r,k max ) = N k (r,g 7 (k max )), A 7 (/c) = A^fc) for 7 ^ 1 and Ai(fc) = 
(o>k/ a oo)Ai(k). In Eq. (9), the role of time is played by the variable /c max . 
The factor g^ in the last term on the right hand side of (9) determines for 
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< 7 < 1 the acceleration of introducing new nodes (of degree one), as /c max 
increases (the number of new nodes introduced in successive equal intervals of 
fc max is for < 7 < 1 a growing function of fc max ). Note that 5^,(&; max )=l for 
7 > 1. The rate equation (9) can be solved using the substitutions (cf. Ref. 20) 
N k (r, /c max ) = rcfc0 7 (fcmax) for < 7 < 1 and N k (l] k max ) = Jfe^ _(A;_1)7 (A; max ) 
for 7 > 1. Thus, according to (7), the systems of equations (5) and (9) are 
equivalent and, consequently, solutions to them yield the same degree distri- 
bution. However, instead of the original attachment rate, determined by Eqs. 
(l)-(3), the rate occurring in Eq. (9) is of the form (4). These rates satisfy 
for 7 < 1 the inequality n 7 (k) < iI 7 (A;), k = 1, 2, except for cases of some 
a>k < floo (when 7 = 1), for which the inequality holds only for sufficiently 
large k. Therefore, using the preferential rate (4) in numerical simulations of 
evolving binary networks enables one to speed up computations, with regard 
to analogous computations based on employing the rate (1), for the same ex- 
ponent 7 < 1 and for the same set of amplitudes {a k } (when 7 = 1). The 
speed up ratio of computations at a given stage of the network generation 
(i.e. at a given moment t) is determined by the quotient of the attachment 
probabilities (4) and (1), or by the quotient of M 7 (t) and /Cmax(7> By means 
of (8), this momentary ratio is given by f 7 (t) = g' 1 (k inax ('j;t)). Clearly, f 7 (i) 
is equal to the evolution acceleration factor occurring in the last term on the 
right hand side of Eq. (9). The total speeding up ratio T 7 (iV), determining the 
reduction of the entire period of time during which a network is numerically 
generated, beginning with one node at the initial time ti — 1 and ending with 
N nodes at the final time tf = N, can be estimated using the relation 



N 



r 7 (iV) = (N-l)- 1 ! g^k max (r,t))dt . (10) 
1 

By applying (7), one finds for < 7 < 1 and for large N that 

r 7 (iV)~0 7 (fc max ( 7 ;i\O) , (11) 

and hence 

f (\nN)~^N , < 7 < 1 , 
N~ , 7 = 1 . 



Notice that, in the case of 7 > 1, Eqs. (10) and (7) imply that r 7 (iV) = 
1. Because r 7 is an increasing function of iV for each 7 from the interval 
[0,1], the application of the rate (4) in numerical calculations in cases of 
< 7 < 1 becomes more and more effective as iV grows (compared with 
analogous calculations with the use of the rate (1)). Clearly, this is of rather 



5 



little importance, as Eq. (5) are exactly solvable for all k when 7 = 1 and 
for asymptotically large k when < 7 < 1. It turns out, however, that, 
in simulations of some networks with more complex normalized linking rules 
than given by Eq. (1), the replacement of normalized attaching rate with 
an unnormalized rate of the type (4) is also possible [22]. Furthermore, the 
application of the linking rate (4) can be very useful for numerical simulations 
of weighted networks with connectivity-driven topologies, especially in cases 
when the weight distribution cannot be found exactly. 

Time needed to simulate networks evolved by node attaching with the proba- 
bility rate with the strictly linear kernel can also be reduced by using, instead 
of the attachment rate (1), a uniform initial attachment rate (correspond- 
ing to 7 = 0), combined with the redirection, with probability r = 0.5, of 
newly created links to ancestors of selected preexisting nodes [21]. Another 
way to enhance the efficiency of simulations of growing networks is the appli- 
cation of parallel algorithms [23]. Although by using different evolution rules 
one can generate in the large N limit networks revealing the same (at least 
for large connectivities) degree distribution, these networks can exhibit dif- 
ferent histories of growth. However, the statistics of weights introduced in 
such connectivity-driven networks cannot necessarily be independent of his- 
tories of the network growth. This can happen when loads are ascribed to 
links as soon as the links are created, and when loads magnitudes depend 
strongly on details of network histories. Consequently, weighted networks can 
display different weight distributions, even if they undergo the same rule of 
link loading and even if they exhibit the same degree distribution. Therefore, 
the rate equation approach to study the weight distribution in networks with 
connectivity- governed topologies should, in general, incorporate rate equations 
for the evolution of the degree distribution. 



3 Evolution of weighted networks 

The process of link loading can be carried out in many different ways, accord- 
ing to various functions performed by real networks [15-19]. For example, to 
mimic functional properties of transportation or communication networks, an 
evolving model has been elaborated by assuming that the topology of the un- 
derlying skeleton network is entirely determined by total weights concentrated 
at vertices to which new vertices are connected, and by assuming that the pro- 
cedure of ascribing a weight to a newly created link involves a redistribution 
(perturbation) of weights, assigned to all links concentrated at the existing 
vertex to which new vertex has been connected [15-17]. However, there ex- 
ist also real networks whose topologies are governed by vertex degrees only. 
Perhaps the most representative example of such networks is the WWW. The 
appearance of a new link in this network depends on the popularity of a given 
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node (webpage) that connects a new node rather than depends on the node 
weight (the total size of documents taken from the preexisting webpage). Nat- 
urally, the popularity of a given webpage can be measured by its connectivity 
(the number of directed hyperlinks joining this webpage with other webpages). 
Although traffic jams in the WWW reduce the speed of sending documents 
from one webpage to others, the total size (usually unknown) of documents 
sent from a given webpage has no effect on the mere act of creating a hyperlink 
between this webpage and another one. Thus, the topology of the WWW can 
be considered as independent of weights assigned to links. 

Here, weighted networks in which node attaching is determined by the prob- 
ability rate (1) (or by the rate (4)) are studied. For simplicity, only the lin- 
ear case 7 = 1, a k = 1, k — 1,2,..., of the preference rate is discussed. 
Each of the links i<->j of a given network is taken to be loaded by a weight 
Wij > 0, and, consequently, each vertex is assumed to concentrate a total 
weight Wi = J2j w i,j- The number of nodes with k links and the total weight 
w is given by 

N k (t) 

N k (w,t)= ]T S(w- Wi ) , (13) 
i=i 



where the sum runs over all N k (t) nodes with k links and 8 denotes the Dirac 
delta function. The rate equations for the evolution of the number N k (w, t) of 
nodes with k links and with the total weight w can be expressed in the form 



d_ 
dt 



N k (w,t) = 

k-i r~ Ek " k 

Pk-i(x) A^_i(tu - x,t) dx - — — N k (w,t) + 



M(t)J ; K " 1V ' ' M(t) 



Sk-l 



5 kyl p k (w - E k + e k ) , (14) 

where M(t) = Mi(t) = ^2jjNj(t), e k ^i is the minimal weight that can be 
introduced to any vertex with k — 1 links in consequence of connecting a new 
link, E k _i denotes a minimal of total weights concentrated at vertices of the 
degree k — 1, and p k (x) is the probability distribution of weights that can be 
assigned to a new link connected to any of vertices of the degree k. Similarly 
to the limit weights e k and E k , the distribution p k (x) is determined by rules 
which govern the weight loading [19]. The first term on the right hand side 
of Eq. (14) is responsible for different processes of the increase of total node 
weights from the value w — x to w, as a result of connecting new links with 
weights x G [e k ^.i,w — E k -i\ (note that w > E k _i + e k ^\ = E k ). The second 
and third terms in (14) are the loss and the new site terms, respectively [20]. 
Obviously, integrating Eq. (14) over all w one recovers the rate equation (5) 
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for the connectivity distribution, in the case of the linear connection kernel. 
In accordance with the argumentation presented in the previous section, the 
evolution of N k (w,t) can alternatively be investigated with respect to k max . 

Since M(t) = 2t [20], one can substitute N k (w,t) = tn k (w). Hence, using (14), 
one obtains 



nk<yW ^ = ~k~T^ I Pk-i( x ) n k-i(w -x)dx + 



w-E k _ 1 



£fc-l 



k + 2 



5 k ,ipk(w - E k + e k ) . 



(15) 



To show that there exists for large w 3> e k ,E k a solution to Eq. (15) of the 
scaled form n k {w) ~ w~ a , a > 1, assume that 

n k (w) ~ n k c k (a) w~ a , (16) 



where 



Uk ~ k(k + l)(k + 2) (17) 



is the exact solution of the rate equations for the connectivity distribution in 
the case of the linear connection kernel [20], and c k (a) is the normalization 
constant, determined by the condition 

oo 

c k (a) J x~ a dx = 1 . (18) 



The scaling expression (16) refers to cases when a weight w is assigned to each 
link with probability determined for all k by the same tail contribution p(w) 
to the distribution p k (w), i.e. p k (w) ~ p(w) forw^> 1, where 

p(w) = \w- a (19) 



with A > being independent of w and equal for all k. Consequently, substi- 
tuting (16) into (15), using Eq. (17) and employing the condition (18), one 
finds for asymptotically large w that 

( i+ £r- i+ (£r ,+ o- ^ 
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where the relation Ek = Ek-i + Ek-i has been applied. Eq. (20) holds if a = 2 
(for w — > oo). Thus, in the case of a = 2, the scaling function (16) indeed 
satisfies Eq. (15) when w tends to infinity. Accordingly, the total weight distri- 
bution P(w) = J2k n k{ w ) undergoes for asymptotically large w the power- law 

P(w) ~ bw~ 2 , (21) 

where the constant b = J2k n kCk{&)- By means of (17) and (18), this constant 
is finite when k~ 2 Ek — > as k — > oo. The weight distribution has also been 
shown to display a power-law behavior in cases of networks in which the growth 
is governed by homogeneous and inhomogeneous attachment rules depending 
linearly on the degree of the existing node that connects new node, and in 
which links receive weights being powers of products of degrees of joined nodes 
[24]. 

The above analysis shows that, in networks with connectivity-driven topology 
and with a special distribution Pk(x) of weights assigned to links, the rate 
equations for the time evolution of the distribution of weights concentrated 
at vertices with a given connectivity have, for asymptotically large weights, 
a solution of the scale-free form. This solution refers to a simple case when 
the numbers of vertices of different degrees are all subjected to the power-law 
for the weight with the same scaling exponent a. Thus, the solution to the 
considered rate equations possesses the power-law form with the same scaling 
index also in cases when the preference rate of attaching vertices is not linear 
[20], provided that the distribution of weights assigned to links is of the form 
(19). The question of whether the studied rate equations for the total weight 
distribution have also other solutions (in cases of other distributions Pk{x) 
than the one considered above), for which the total weight distribution P(w) 
would reveal a scaling form at large w, is open. Clearly, what values could 
take the scaling index describing the power-law decay of the total weight 
distribution, or whether this distribution could undergo the power-law at all, 
it would depend on specific rules of link loading in an individual network. Note 
that the rate equation approach has been applied to investigate weight and 
degree distributions in networks with weight-driven topologies [25]. 
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